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Abstract 

The process of pair creation by a photon in a strong magnetic field is 
investigated basing on the polarization operator in the field. The total prob- 
ability of the process is found in a relatively simple form. The probability 
exhibits a "saw-tooth" pattern because of divergences arising when the elec- 
tron and positron are created at threshold of the Landau energy levels. The 
pattern will be washed out at averaging over any smooth photon energy dis- 
tribution. The new results are obtained in the scope of the quasiclassical 
approach: 1) in the case when the magnetic field B <C Bq, (Bq is the critical 
field) the new formulation extends the photon energy interval to the case 
when the created particles are not ultrarelativistic; 2) the correction to the 
standard quasiclassical approximation is found showing the range of appli- 
cability of the approach at high photon energy as well. The very important 
conclusion is that for both cases BcBo and B > Bq the results of the qua- 
siclassical calculation are very close to averaged probabilities of exact theory 
in a very wide range of photon energies. The quasiclassical approximation is 
valid also for the energy distribution if the electron and positron are created 
on enough high levels. 
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1 Introduction 



Study of QED processes in a strong external field is stimulated essentially by exis- 
tence of very strong magnetic fields in nature. There is a series of evidences for the 
existence of neutron stars possessing magnetic fields close the critical field strength 
B = m 2 /e = (m 2 c 3 /eh) = 4.41 • 10 13 G. The rotating magnetic dipole model, in 
which the pulsar loses rotational energy through the magnetic dipole radiation, was 
confirmed with the discovery that the spin-down power predicted for many pulsars 
is in a quite good energetic agreement with the observed radiation, and gives mag- 
netic fields B ~ 10 11 — 10 13 G [lj. There are presently around 1600 spin-powered 
radio pulsars. 

Another class of neutron stars was discovered at X-ray and 7-ray energies and 
may possess even stronger surface magnetic fields B ~ 10 14 — 10 15 G. Such stars are 
now referred to as magnetars [2j, since they most probably derive their power from 
their magnetic fields rather than from spin-down energy loss. Two different classes 
of objects are thought to be magnetars: the Soft Gamma-ray Repeaters (SGR), dis- 
covered as sources of short bursts of hard X-rays (> 30 keV) with super-Eddington 
luminosity and the Anomalous X-ray Pulsars (APS), discovered as persistent, soft 
(< 10 keV) X- ray sources with pulsations of several seconds and spinning-down 
on time scales of 10 4 — 10 5 years. An increasing number of common properties, 
pointing to a close relationship of these two apparently different classes of objects 
and leading to their interpretation as magnetars, has been found (see e.g. [3]). 

Pair creation by a photon in a strong magnetic field is the basic QED reaction 
which can play the significant role in processes in vicinity of a neutron star. 

Investigation of pair creation by a photon in a strong magnetic field was started 
in 1952 independently by Klepikov and Toll [HE]. In Klepikov's paper [6], which 
was based on the solution of the Dirac equation in a constant and homogeneous 
magnetic field, the probabilities of radiation from an electron and e~e + pair creation 
by a photon were obtained for magnetic field of arbitrary strength on the mass 
shelfl (k 2 = 0, k is the 4-momentum of photon). In 1971 Adler [7] calculated the 
photon polarization operator in the mentioned magnetic field using the proper- 
time technique developed by Schwinger [BJ and Batalin and Shabad [9] calculated 
the photon polarization operator in a constant and homogeneous electromagnetic 
field for k 2 7^ using the Green function in this field found by Schwinger [8]. In 
1975 Strakhovenko and present authors calculated the contribution of a charged- 
particles loop with n external photon lines having applied the proper-time method 
in a constant and homogeneous electromagnetic field [ID]. The explicit expression 
for the photon polarization operator was found. In a purely magnetic field (E=0) 
the imaginary part of the polarization operator, which determines the probability 
of the pair creation in general case (k 2 7^ 0), was analyzed in detail near the 
root singularities originating at each threshold of e~e + pair creation in the state 
with the given Landau quantum numbers. This peculiarity of the pair creation 

1 We use the system of units with Ti — c — 1 and the metric ab = a^b^ = a°b Q — ab 
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process was observed by Klepikov [6] . Shabad [TT] performed an extensive analysis 
of the photon polarization operator found in [5] in a purely magnetic field. It 
was shown that its imaginary part in the limit k 2 — > reduces to the result of 
Klepikov's calculation [6]. Daugherty and Harding [12J analyzed in detail the total 
probabilities (the attenuation coefficients) of the pair creation by a polarized photon 
in the superstrong field B ~ Bq basing on the calculations of Toll and Klepikov and 
paying the special attention to the threshold region u > 2m, where the probabilities 
exhibit a "saw-tooth" pattern due to the mentioned root singularities. 

For high Landau levels Klepikov [6] (and Toll) found the quasiclassical repre- 
sentations for the probabilities of radiation and pair creation in the case B <C Bq. 
In 1968 the present authors developed the operator quasiclassical method for the 
consideration of electromagnetic processes in an external field [T3] which is valid 
also in nonuniform and (or) depending on time field under conditions: 1) B <C B Q , 
2) the charged particles have relativistic energies. Both the radiation and pair 
creation by a photon were considered. The method is given in detail in [IH [15]. 
Tsai and Erber [16J deduced the quasiclassical probability of pair creation using 
the photon polarization operator in a magnetic field calculated by Adler [TJ. 

In the Sec. 2 the exact probability of the pair creation by a photon was ob- 
tained for general case k 2 ^ starting from the polarization operator in a mag- 
netic field. The technical details of the derivation are given in Appendix A. In 
Sec. 3 the standard quasiclassical approximation (SQA) is outlined in the limit 
H = B/Bq <C 1, r = oj 2 /Am 2 3> 1, where uj is the photon energy. The quasi- 
classical characteristics of the process depend on the parameter k = 2/i-y/r. The 
correction to SQA is found. Details of calculation of the correction starting from 
the polarization operator in a magnetic field is given in Appendix C. If the photon 
energy is not high, so that /i <C r — 1 <C /i -2 and SQA is not applicable, the new qua- 
siclassical approach (at low energy near threshold (th)) is developed in Appendix 
B. One can expect that the photon energy distribution in a vicinity of neutron star 
is wide and smooth. At averaging over this distribution the "saw-tooth" pattern 
will be washed out. It is very important that the result the averaging is very close 
to SQA even at /x > 1. The corresponding analysis has been carried out in Sec. 3. 
The energy spectrum of created pair is discussed in Sec. 4. Here SQA is very useful 
also. A remark concerning derivation of the energy spectrum in SQA is given in 
Appendix D. 

2 Probability of pair creation by a photon: exact 
theory 

Our analysis is based on the expression for the polarization operator obtained 
in [ID], see Eqs.(3.19), (3.33). For pure magnetic field (E=0) this polarization 
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operation can be written in the diagonal form 

k^k y 

rp" = -$>/W, hh -<\ r ptk = o, ^^ = — -g^, (2.1) 



where 



k ±y /k?(u 2 - hi) 1 Bk ± " B^u 2 - k\ 
«i = f)i(r — g), k 2 = Ki - qfl 2 , K 3 — Ki+rQ s , 

r _ " 2 - fcj g= iL r _ g = _^ (22) 
4m 2 ' 4m 2 ' 4m 2 ' 

where the axis 3 directed along the magnetic field B, k^B = 0, k± = \J\<?\_ 1 oj is the 
photon energy, F^ u is the tensor of electromagnetic field, F*^ v is the dual tensor 
of electromagnetic field, 

9 1 00— iO 

corn f f 

Qi = dv / fi(v,x) exp(iip(v,x))dx. (2.3) 



where 



7T 

-1 

Here 

cos vx cos x sin vx cos v x — cos x 

fi[v,x) = — v — , f 2 {v,x) = 2 — fi{v,x), 

sm x sin x sin x 

p i \ / \ n / \ 1 / \ COS Ob COS VOC 7 / \ / \ 

J3V0, x) — (1 — v ) cot x — Ji(v, x), ip(v, x) = z bivXxvzA) 

smx 

2g N 1 -r(l -v 2 ) B . , 

Z = JL, h {v) = ^ L, fi = —. 2.5 

fj, [i B 

Let us note that the positive x axis in the integral Eq. fl2.3p is turned slightly down, 
and in the function Qi in the integral over x the subtraction at \i — is implied. 

The real part of the polarization operator defines the dispersive properties of 
the space region with magnetic field. At r < 1 {b > 0) one can turn the integration 
contour over x to the negative imaginary axis (x — > —ix) so that 

am 2 r °f . . / cosh x — cosh vx\ 

Hi = dv ifdv, —ix) exp — ox — z dx. (2.6) 

n J J \ srnhx / 

-10 v 7 

The functions 0« in Eq. (12.61) are real. These functions, studied in [T7] at u < 2m 
and arbitrary value of //, define the index of refraction 712,3. On mass shell k 2 = 
one has 

<s = 1 - ^- (2-7) 
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For the mode 2 the photon polarization vector e 2 is perpendicular to the plane 
formed by the vectors B and k and for mode 3 the photon polarization vector e 3 
lies in this plane (n 2 — > nil and n 3 — > n± in Adler's notation [7]). 

At r > 1 the polarization operator eigenvalue K{ acquires an imaginary part, 
which determines the probability of e~e + pair creation per unit length W% by a 
photon with a given polarization 

W = - Ime^* U^. (2.8) 

Details of calculation basing on Eqs. (12. 21) . (12.31) are given in Appendix A. Using 
e i — Pi we get the explicit expressions for Wi {i = 1, 2, 3) at k 2 ^ 0: 

Wi = = E f 1 - —) m^rA (2-9) 

where 5„o is the Kronecker symbol, i?(gf) is the Heaviside function: "d{g) = 1 for 
9>0, (<?) = for <? < 0, 

m + n m — n n 2 /i 2 

/ = ^— > /c = ^— > = r-l-m// + — , 

U = Kr)], d(r)= 2(r ~^ , (2.10) 

/i 

Here [d(r)] is the integer part of d(r), 

d i=\Z- X ) G ' rf 2 = -G + 4/iW(A; - l)V^\{z)L^\z), 
d 3 = (l + ^ - ^ F + 2^(* - l^L^), 
C = {^f - «£pj F - 2»W(k - WtiW, 

F=(^(z)) 2 + ^-l)i-(LL 1 (^)) 2 , (2.11) 

where L^(^) is the generalized Laguerre polynomial, 2 is defined Eq. fl2.5p . Obtained 
probabilities Wi agree with found by Shabad [TT] but presented in much more 
compact form. 

On the mass shell (k 2 = 0) in Eqs.( l2.9l) - (l2.11j) one has to put r = q. In this 
case W\ = and only two polarizations i = 2 and i = 3 remain. The probability 
of pair creation averaged over the photon polarizations acquires especially simple 
form (in the frame where ^3 = 0) 



2 A c 2 ; y/gl\ 



[2.12] 



The probabilities of pair creation by a photon Eqs. (12.91) . (12.121) contain the 
factor 1/y 7 ^. The function g = 0, if 



where e(l) is the energy of charged particle on the Landau energy level. So the 
probability diverges when the electron and positron are created on the Landau 
levels with electron and positron momentum p% = 0, see e.g. [6l [TO]. The origin 
the singularity is due the properties of the space volume in the lowest order of 
perturbation theory (infinitesimally narrow level). The reason why this singularity 
is not a pole but a branch point is that motion along a field is not quantized. 

3 Probability of pair creation by a photon: 
quasiclassical approximation 

In the case /z-Cl, r — 1 ^> 1 (l ma x 3> 1) the mentioned above divergences are 
situated very often: A/ ~ Ar/fi ~ 1. On the photon energy scale the characteristic 
distance between the mentioned peaks will be Au ~ ujq, where uo is the frequency 
of motion of created particle in a magnetic field. This frequency defines the dis- 
tance between energy levels in the quasiclassical approximation. So the resonance 
properties of pair creation probability become apparent if the width of energy lev- 
els and the effective spread of photon energy are small comparing with ujq value. 
In the opposite case the spectrum peculiarities will be washed out and the process 
probability will be quite adequate described by approximate expressions ultimately 
connected with quasiclassical character of created particle motion at large Landau 
quantum numbers. 

When the magnetic field is weak comparing with critical one (/i < 1) and 
for high energy photons (r > 1) the quasiclassical approximation was created in 
[HJ [13], [16], where the process probabilities were derived using the different ap- 
proaches. In [B] the asymptotic of Eq. fl2.12p were obtained. The approach in [T3] 
was developed for jj, <C 1 and the particle motion was considered to be quasiclassi- 
cal (and relativistic) from the very beginning; it is applicable for nonuniform (and 
non-stationary) fields. The pair creation probabilities in the form Eq.( lC.6l) (see 
Appendix C) were found in [16] . Below we will call it the standard quasiclassical 
approximation (SQA). 

In Appendix C the corrections to SQA were calculated. For probability of pair 
creation by unpolarized photon the correction is 




(2.13) 



2 2 

am jj, 



i 



dv 



60^ 



I 



2(1 + v 2 - 27z 2 )K 1/3 (z) + 3(7 - v 2 )zK 2/3 {z) 



1 — v 



2 ' 



(3.1) 



6 



where 

z = Tr\ 2W K = ( 3 ' 2 ) 

3(1 — V 2 )K 

At k <C 1 the main contribution is given by the term oc z 2 , it is 

^ " Book 2 ex P(, 3 J> w " i5K 3- V- 6 ) 

It is seen that at k <C 1 the SQA is applicable if // 2 /k 3 <C 1. 

At r < /i -2 / 3 SQA becomes inapplicable. In this region one can use the process 
probability found in Appendix B, which is valid if // <C r— 1 <C //~ 2 . Using Eq.( lB.5l) 
we get for the probability of pair creation by unpolarized photon 

W™ = , 3r ~ 1 exp , (3.4) 

4w yJr(r-l)l(r)P(r) \ V J 



l(r) = ln V ^ + ? , /3(r) = 2^-(r-l)Z(r). (3.5) 
/r — 1 



where 



For large r> 1, limited by the condition r <C // 2 (k < 1), we have 

At K 3 ^> // 2 Eq. (13.61) agrees with Eq. (|3.3|) . The last expression has essentially wider 
region of applicability than pure W^ S ® A '. 

The behavior of pair creation probability for unpolarized photon at fi = 0.1(1/// ^> 
1) as a function of r is illustrated in Fig.l. The probability in saw-tooth form is 
calculated according to Eq. ( 12. 12ft . It is seen that the characteristic distance be- 
tween peaks is of the order // = 0.1. The result of averaging over this interval 
(— ji/2 + r-=-///2 + r) is given by the dashed line which is smoothed out all pecu- 
liarities of original behavior and which is in a very good agreement with the thick 
curve calculated according to Eq. (13.41) . 

At k > 1 Eq. (13.41) becomes inapplicable. However in this case SQA is valid. 
From Eq. (1C.6P we have the probability of pair creation by unpolarized photon 

w <*A> = J!g_f«Z± K 2/3(z)dv , (3.7) 

where z is defined in Eq. ( 13.21) . The correction to this probability is given by 
Eq. ()3.ip . its relative value at k > 1 is < // 2 . The ratio W^ tfl) and W( SQA) is shown 
in Fig. 2. The curves 2 and 3 are for polarized photons, the thick curves is for 
unpolarized photon (see Eq. (13.41) and Eq. (j3.70 ). The point where the ratio attains 
1 is the boundary of the applicability region of the corresponding approximation. 
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The relative contribution of the correction Eq. fl3.ip to SQA at [i = 3 is shown 
in Fig. 3. The upper curve is the probability of pair creation by a photon in SQA 
Eq. fl3.7p . The lower curve is the sum W^®^ + W^ 1 ' (the correction is negative at 
low r). It is seen that the correction is small in the region of SQA applicability 
r > fi/2 = 1.5. 

It is very important that the quasiclassical approximation appears to be ap- 
plicable also for super-strong magnetic field [i > 1. This is true if the character- 
istic parameter l max ^> 1. It is assumed in this case that parameter k is large 
(k = 2^fr[i 3> 1). In the limit k ^> 1 the asymptotic expansion of the probability 
Eq.fJSZD is 

(mA , am 2 ( 2tt \ 3 7 / 6 5 r 3 (2/3) . , 

w is W = _^ - T + ...), B = __i2_J = 1 . 1925 .. (3.8) 

The asymptotic expansion of correction Eq. fl3.ip is 



WV = —^^(AK^-6n + ...), A= " K £ '*X'Z =8.191.. (3.9) 



2 



? f,.Ms c _, \ A 3^213(1/3) 



iroo QOVSk v " ' ' 5 T(2/3) 

The relative value of the correction defines the boundary of region where SQA is 
valid for the super-strong magnetic field [i > 1 

W^L = 3-r/s _7_T^!3)jl_ L \ 066^^1-^ 

^(SQA) ° 125 r 4( 2 /3) K 4/3 y 1 Ak 1/3 • • ^ " U - UU ° K 4/3 ^ R l/3 / 

(3.10) 

It is seen from this expression that at [l > 1 the applicability of quasiclassical 
approximation is controlled by the parameter (/i/2r) 2 / 3 = /~^ 3 1. So, both 
in the case ^ <C 1 and in the case /i > 1 the quasiclassical approximation is 
valid when the particles of pair are created in states with large Landau quantum 
numbers. It should be noted concerning Eq. f)3.10p the smallness of numerical factor 
and relatively large contribution of the second term oc k^ 1 ^ 3 . This means that the 
correction is small even at low l ma x, and one have to use the exact expression 
Eq. fl3.ip . Under these circumstances one can use Eq. (13.71) with the correction 
Eq. fl3.ip starting from r ~ /x. This is illustrated in Fig.4 where the probability of 
pair creation is shown at fi = 3. The probability in saw-tooth form is calculated 
according to Eq. (12. 12j) . It is seen that the characteristic distance between peaks is 
of the order fi (here fi = 3) just as in the case shown in Fig.l where the probability 
for [i = Q.l was presented. The result of averaging over this interval {—fi/2 + 
r 4- fi/2 + r) is given by the dashed line which is smoothed out all peculiarities of 
original behavior and which is in a very good agreement with the thick curve which 
is quasiclassical approximation calculated according to Eqs. (13.71) and (13.11) . 



4 Energy distribution of created pair 

The energy distribution of created pair is symmetric with respect to exchange of 
electron and positron, i.e. it is symmetric respect energy uj/2 and one can consider 
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this distribution in interval x = e^/uj > 1/2. In the frame where k 3 = the energy 
of electron and positron is 



MO = y/li + m 2 {l + 2fjl), e+{k) = ^Jp 2 3 + m 2 {l + 2fik), u = £_(/)+ e+(k). 

(4-1) 

From Eq. fl4.ll) follows that p 2 /m 2 = g, (g is defined in Eq. fl2.10p ). In the plane (k,l) 
the allowed states of electron and positron are bounded by the axes and the curve 
pi = which is parabola with the nearest point to origin l n = k n = [(r — l)/2/x]. 
When k — the maximal value I = l max (l ma x is defined in Eq.( 12.10j) ). Substituting 
p 2 /m 2 in £_ we get 

x = — = -(l + ^£), n = l-k. (4.2) 



u 2 V 2r 

So the distribution over energy of created particle is of the discrete character. The 
probability W(n) is defined by Eq. fl2.12p . where one has to perform summation at 
fixed n = I — k. The expression for spectrum contains the factor 1/ yfg discussed 
above (see Eq. fl2~T3l l 

The energy distribution over electron energy in SQA (see, e.g. [T3], Eq.(3.50)) 

is 



dW^ SQA ^ am 2 



dx y/3lT(jJ 

where 



x 2 + (I - x) 2 



K 2/3 (0 + J K 1/3 (y)dy 



x(l — X) J£ 



(4.3) 



^ 3kx(1 — x) ^ 



At low values of the parameters [i and k the energy distribution of pair creation 
probability has a sharp peak when e_ = e + = u/2 (for details see Sec. 3. 2 in [15J). 
When /i C 1 and k > 1 the SQA is not valid on an edge of the energy spectrum 
where the created particle is no more relativistic (1 — x ~ m/u). However at 
smaller energy (1— x<C 1/k = m/(fiuj)) the differential probability is suppressed 
exponentially, so the contribution of the interval 1 — x ~ m/uj, where SQA is not 
valid, is negligible small. 

For the super-strong field /i > 1 in the region of applicability of the quasiclassical 
approximation (r ^> /i) the parameter k = 2/zi/r ^> 1. In this case the upper limit 
of electron energy xj, — 1 — 1/ (2y / r) is lower (or of the same order) than the electron 
energy in the maximum of the energy distribution of Eq. fl4.3p (x m ~ 1 — 1.6/ k). 
The region of applicability of the quasiclassical approximation near x^ is defined by 
the quasiclassical character of positron (electron) transverse motion (p 3 = 0, k ^> 1) 



\2 x ^ A* i An 



{l-x) 2 - — l-x > An> v ^, (4.5) 



where n m = l max Eq.(l2.1(jp. 

In Fig. 5 the spectrum of created particles is shown for ji = 3, r = 110 {n m = 66). 
The spectrum in saw-tooth form is calculated according to Eq.f l2.12l) . The dashed 
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curve is the spectrum averaged over the interval (— 2/i + r-=-2/i + r). The thick curve 
calculated using Eq. (14.31) . It is seen that this curve is in a quite good agreement 
with the result of averaging. 

In this case the argument £ <C 1 in Eq.f l4.3p and one can perform the expansion 
(k 3 = 0) 



dW (SQA) 



dx 



amy 



Ci 



x 2 + (l 



KX[1 



- xf J_ 

x)y/z + 3k 



Co 



1 



[KX(1 - x)) 5 /3 



+ . . . 



(4.6) 



where 



3V6 
2vr 



r(2/3) = 0.2588 



Co 



3-1/6 

47T 



r(l/3) = 0.1775.. 



(4.7) 



In Fig. 6 the energy distribution of the created electron(positron) at fi = 3 is shown 
for different r: r = 30, n m = 16 (the upper curve), r = 110, n m = 66 (the 
middle curve), r = 500, n m = 318 (the lower curve), calculated according with 
Eq. fl4.3p (the solid curve) and according with Eq. fl4.6p (the dashed curve). Very 
good agreement of corresponding curves is seen in the region of SQA applicability 
{x>n-V\{l-x)>n^). 



5 Conclusion 

In this paper we study the process of pair creation by a photon in a strong magnetic 
field basing on the polarization operator in the field calculated by different methods 
in a set of papers [3, [9j [10]. By using the decomposition of functions into series 
containing the Bessel function (see Eq.( 1A.4p ) we deduced the imaginary part of the 
polarization operator (which is the process probability) into the relatively simple 
combination of the generalized Laguerre polynomial (see Eq. ()2.12j) ). This proba- 
bility exhibits a "saw-tooth" pattern because of the factor 1/ \fg. This pattern will 
be washed out at averaging over any smooth photon energy distribution. 

The range of applicability of the quasiclassical approximation extended in the 
photon energy interval where the created particles are no more ultrarelativistic (see 
Appendix B). Found correction to standard quasiclassical approximation (SQA) 
permitted to extend the range of applicability of the quasiclassical approximation 
to the region r ~ // at // > 1. From the performed analysis follows the remarkable 
conclusion: for both < 1 and ji > 1 the results of the quasiclassical calculation 
are very close to averaged probabilities of exact theory in very wide interval of pho- 
ton energies. For the total probability of pair creation the quasiclassical approach 
is valid if (see Eq.fjnOD) 

l max = [ 2(r ~^ l » 1 (5.8) 
H 
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In relatively weak field (/i <C 1) it is valid not far from threshold: (r — 3> 1. 
For this energy interval the relative correction to the probability Eq. (I3.4|) is of the 
order /~^, 2 . In superstrong field /i > 1, for relativistic motion of the particles of the 
created pair and if l max ~ 2r/fi ^> 1 the SQA is valid. The expansion parameter 
in this case is I^ 3 (see Eq. (IC.llT) ). 

For such field and energy it is helpful to compare the probability of pair creation 
by a photon W with the radiation length L ra( j for photon emission process found 
in [IS]. Conserving the main terms of the decomposition over the Landau energy 
number / (see Eq. (12.13]) ) in the SQA as well as in the first correction we get (k 3 = 0) 



r-i 

^rad 



11/6 + 



w = a -f 

Li 



1 



+ 



,1/6 1 ,5/6 
I'liiax lmax 



(5.9) 



where 



32 



B\ 1/2 

e 



a 



e 

foe 



1 

137' 



18)" 1/6 r(2/3) = 0.33045... 



(12)' 1 / 3 r(l/3) 
20 T(2/3) 



0.043206..., b r , 



15 n rl/6 r 2 (2/3) 
' 7 1 ; r 2 (i/3) 

7(12)-i/3 r 4 (l/3) 



0.33819..., 



1125 T 4 (2/3) 



0.041633..., (5.10) 



where the numbers are taken from Eq.(46) in [18] and Eqs.( 13.8l l). (13. 9jl ) above. 
It should be noted that the physics characteristics in Eq. (15.91) ) don't depend on 
the electron mass. It's important that in the formula not only coefficients a 7 
and a p are very close, but the same properties possess the coefficients 6 7 and b p 
reflecting proximity of details of mechanisms of radiation and pair creation by a 
photon in a magnetic field. One can expect that this property will be essential for 
consideration of electron-photon shower in a magnetic field. Set of kinetic equations 
for such shower with the particle polarization taken into account was investigated 
in [T9l[20]. 
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A Appendix 

Let us consider the integral in Eq. fl2.3j) 



1 oo— iO 

%= Idv I /,(,,«) «p(^(„,,))d,. (A.l) 



-1 



Using the formula 



- cot x^-e^ = ( f 2 - f x - — cot x] e^, (A.2) 
iz ax \ z I 



and integrating by parts the integral over x in Eq.( 1A.ll) we get 

oo— iO oo— iO oo—iO 

f f 2 e^dx= [ f 1 e i *dx + - I (bcotx L-) e ^dx (A.3) 

J J z J V sin xJ 



Using Eq. 7.2.4. (27) [21J we have 
e" rfcosra = E ( 2 " <W (-i)"^nW cos(n ra ), * 



„=o sina; ' 

e -i tc osvx CQSVX = ^-ttcosvx = J2 \ l - ^.j (J n+1 (t) - J n ^{t))c0s{nvx), 

Id 2? 00 

e - it cosvx sinvx = -ttcosvx = _ y (_ ? )« n J n (t) sin ( nra ) 5 (A.4) 

itv dx t i 

where J n (t) is the Bessel function. Taking into account that the functions fi and 
ip in the integral over v in Eq. (12.31) are the even functions of v, we can perform 
following substitutions in sums in Eq. (lA.4p 



cos(nvx) -> e inra , sin(nra) — > -ze irira . (A. 5) 

So the expression for T{ Eq.f lA.ll) can be written as 



1 oo— iO 



Tj = E^i ) ^ — / ^ / Fi n \ v i x ) exp(ia n (v)x)dx, (A. 6) 

10 

where 



Ff = (-i) n exp(«2cotx) 



1 ~ "7T J (*) - - — COt Xj n (t) 



SIM \ 2 / 2 



F (n) = F (n) + eX p(^C0tx) 2(2 ^ (fecotz - J n (t), 

z V sin x / 

F 3 (n) = (-i) n exp(^cotx)(l - v 2 ) cotx (2 - <5 n0 ) J n (t) - F[ n \ 
a n (v) = nv — b, (A. 7) 
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(n) 

and the function F{ is a periodic function of x. 

Let us note that at x — > —zoo the asymptotic of J n (z/ sinx) is 



J n — — ~ J n (2tze- X ) ~ ^-e"™". (A.8) 
Vsmnx/ n! 

In the case a n (t>) < n the integration contour over x in Eq. l1A.6l) can be unrolled 
to the lower imaginary semiaxis and the contribution to the imaginary part of 7^ 
vanishes. 

One can represent in the form 

\ 2 r 

T/ n) = dv Ft\v,x)exp{ta n {v)x)^2e 2mkan(v) dx 



-1 o 
1 2tt 

dv 



k=0 



I _ e 2wia n {v) _|_ 



F^ n \v, x) exp(ia n (v )x)dx. (A. 9) 



We will use in Eq.( lA.9l) the known equality 

1 V 



\ g27na n (i>) _)_ ^0 1 f>2wia n (v) 

By virtue of the made above observation 



inS(l~e 27ria " {v) ). (A. 10) 



- m5(l - e 2 ™ n{v) ) = -in ^5(1 _ e ^(M»)-m)) _, _ ^5(a n (u)-m). (A.ll) 

ti m~>n 

Substituting Eq. flA.lOj) and Eq. flA.l If) into Eq. (IA.9[) and taking into account that 

Ft\v,x + n) = (-irFt ) (v,x), (A.12) 



we have 



rp{n) 
i 



1 d 

" 1 ) n o^ > / ttt / i 71 / n) (^,a;)exp(za n (u)x)rfx 

2 J smyKa n {v)) J 



+ 53 51 — 2|a' fa I — J F i n \v h 2,x)exp(imx)dx, (A.13) 



where 



n// m 2 /i 2 1 + m// 

1,2 2r V 4r 2 r 



, / , , 2 r- n 2 /i 2 . . 

a n( u ) = — V r 9\ n i m ; r ); gin, m, r) = r — 1 — mfi -\ . (A. 14) 

/i 4r 
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It is seen from Eq. (1A.13j) that n + m is even. Using Eq. (IA.120 one can represent 
the imaginary part of T^ n ' m ^ in the form 

tt/2 

ImT/ n '"^ = — i = / Fj; n \vi2,x) exp(imx)dx (A.15) 

vi, 2 ^rg(n,m,r)_J 



/2 



Along with integers m and n we will use also I = (m + n)/2 and k — (m — n)/2. 
We turn now to calculation of integrals in Eq. (lA.15p . Let us consider the integral 



C^Jz) = (-l) n+1 [ J n+1 ( —) e i(mx+Zcotx)J^_ (A. 16) 

J \sinxj sinx 

-tt/2 

Changing the variable t = cot x we get 



CiW = H) B+1 / ^n+i (^) W2 eiZt -7^=j- ( A - 17 ) 



One can close the integration contour in Eq. flA. 17|) by an infinite half-circle in the 
upper half-plane and than contract the contour to the point t = i. So, we have to 
calculate the integral in Eq. (1A.17l) over the contour Changing successive the 

variables t — > it, y = (t — l)/(t + 1) we get 

C™ +l {z) = e~' I J n+1 y^e-Vd-ri fy (A.18) 



(o+) 

Using Eq.l0.12.(20) 
1 



y/y{i-y)' 



^^^=£^MM)V, (A.19) 

where I n (x) is the modified Bessel function, L^(z) is the is the generalized Laguerre 
polynomial, we get 

C™ +l {z) = 2me- z z n+1 ^^{Llt{{z))H{k - 1), n = I - k. (A.20) 
In the same manner one can calculate the integral 

7r/2 



S™(z) = 2(-i) n [ cotxJ n (—)e i{mx+zcotx) dx 
J \smxj 



-tt/2 

h\ / 

2irie- z z n jF(l, k, z), F(l, k, z) = (L n k (z)) 2 + -(L^z)) 2 ^ - 1). (A.21) 
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Somewhat different structure has the integral 



tt/2 



2ie 



j j Z ^ c i(mx+z cot x) dx 



SIM 



-tt/2 



4 +1 ( ^ ) f m/2 e-* 



(o+) 



EE 



(o+) 

Using Eq.10. 12.(9) [21]: 



(L n (z)) 2 v j+S ~ k 



sin a; 
dy 

-47ie~ z z 



k-i ■ 



S(JT^ ( ^ W)2(A - 22) 



E 



one obtains 



L-(x)LUy) 



k\ 



(k + n — 1)! x — y 



(A.23) 



-47re~ 2 ^' 



(/-I)! 



Lt{z)j- z Ll_ x {z) - LUiz^L^z) 



k} 

A7ie- z z n -l 



L n k tl(z)Ll-\z)-LUz)LU(z) 



(A.24) 



In derivation of the last expression Eqs. 10. 12. (15), (16), (24) [2lJ was used. In con- 
clusion of this Appendix let us note, that x/ ri ' m ' ) contains the combination 



C^. 1 (z)+C^. 1 (z) = 2ine-'z n - 1 j 



^{Ll+{{z)f + l{Lr\z)T 



2me~ z z n - 
2me- z z n ' 



l\ 



fc 



(A.25) 



It should be noted that the above derivation is the most direct path from the process 
probability in the form of the integral Eqs. fl2.2l) . fl2.3l) to the process probability in 
the form Eq.flZS}. 



B Appendix 

Quasiclassical approximation at low photon energy 

(u ~ m, k 2 = 0, r = q) 

In the field which is weak comparing with the critical field B/B Q = \i <C 1 and 
at relatively low photon energy (r ~ 1) the created particles occupy mainly states 
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with high quantum numbers if the condition r — 1 3> [i is fulfilled. In this case the 
quasiclassical approach is valid, but created particles are no more ultrarelativistic, 
so the standard quasiclassical formulas [131 HU US] are nonapplicable. We will 
develop here another approach, using the method of stationary phase at calculation 
of the imaginary part of the integral over x in Eq. (12.31) . Granting that the large 
parameter 1/fi is the common factor in the phase ip(x) and doesn't contained in 
the equation ip'(x) = which defines the stationary phase point x${r ~ 1) ~ 1. In 
this case the small values of variable v contribute to the integral over v, so that 
one can extend the integration limits to the infinity. So we get 



2 oo oo 



OiXfi P f it r i I 

ImQi ~ i— — I dv I fi(v , x) exp < (p(x) + v 2 x(x) > dx, (B.l) 

— oo — oo I ^ J 

where 

(f(x) = 2rtan ( — ) + (1 — r)x, x( x ) — rx [1 ) . (B.2) 

\2J \ sins/ 

From the equation (p'(x) = we find 

tan(^)=-^=, x (r) = -il(r), /(r) = ln^±| (B.3) 



Substituting these results in the expressions which defines the integral in Eq. (IB. lj) 
we have 



iip(x Q ) = (3(r) = 2y/r - (r - l)Z(r), i<p"(x ) = ix(xa) = —l{r)(3{r) 

if 2 (v = 0,x (r)) = ^i, -if 3 (v = 0,x (r)) = -L (B.4) 



Performing the standard procedure of the method of stationary phase one obtains 
for probability of pair creation by a polarized photon not far from threshold 

In spite of made above assumption r ~ 1 Eq. (IB.5l) is valid also at r ^> 1 if the 
condition j3(r) ^> fi is fulfilled. This can be traced in the derivation of Eq. flB.5l) . 
The first two term of the decomposition of the function j3{r) over power of 1/r are 

Pir) 44 ,„ , 

— — ~ 1 (B 6) 

fji ~ 15/zr 3 / 2 1 ' 

It follows from this formula that that applicability of Eq. (1B.5p is limited by the 
condition r <C yT 2 . If the second term much smaller than unity the exponent with 
it can be expanded. As a result we have from Eq. (1B.5l) at /i~ 2 / 3 < r <C [i~ 2 
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The main term in above expression coincides with the probability of pair creation 
by a photon in standard quasiclassical theory at k = Ifi^fr <C 1. The correction 
in Eq. (lB.7p determines the lower boundary over the photon energy of standard 
approach applicability (k 3 ^> /i 2 ). So exists the overlapping region where both the 
formulated here and the standard approach for high energy are valid. 



C Appendix 

Corrections to the standard quasiclassical approximation 

(k 2 = 0, r = q) 

The standard quasiclassical approximation valid for ultrarelativistic created 
particles (r ^> 1) can be derived from Eqs. fl2.3p . (12.41) by expanding the functions 
/2(f) x), fs(v, x), if)(v, x) over x powers. Taking into account the higher powers of 
x one gets 



f2(v,x) 
fa(v,x) 
ip(v, x) 



1 - v 2 r l 

(3 + v 2 )x + — (15 - 6v 2 - v 4 )x 3 
15 



12 

1-v 2 

r(l — v 2 ) 
12/7 



(3 - v 2 )x + —(15 - 2v 2 + 3v 4 )x 3 
60 



2\2 



30 



x 

x" I . 



(c.i; 



Here the first terms in the brackets give the known probability of the process in the 
standard quasiclassical approximation, while the second terms are the corrections. 
Expanding the term with x 5 in exp(iip(v,x)) and making substitution x = fit one 
finds 

1 00 



Imfli = 

92(v,t) 

where 



.am 2 fx 
2tt 

1-v 2 



12 

1-v 2 



dv j gi(v,t)exp 

1 — 00 

, 9-7 

fit 



t d 



dx. 



(3 + v< 



90 



^H 5 + ^(15 - 6v 2 - v A ) 
15 



6 



-fit 



(3 - v 2 ) - 1 



90 



60 



, \(1-v 2 )k] 2 r- u 

16 m 



(C.2) 
(C.3) 



We will use the known integrals 



OO y g \ 

| cos ft + ^-3] = 



z), z 



-00 

OO 



3v^ 3(1-v 2 )k' 



I tsin ( t+ 4) 



3V3 



z 2 K 2/3 (z) 



(C.4) 
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Differentiating the first integral over £ one gets 

— oo N ' 

<7'-K)--3£te) , < rfr 

4zK 2/3 (z) - + z 2 ) K 1/3 (z 



27^3 3 



z d zK. 



2/3 1 



^1/3 («: 



81^3 

A 

16 



(C.5) 



Substituting the integrals Eqs.( 1C.3j) -( 1C.4l) in Eq.( 1C.2l) we obtain the probabilities 
of pair creation in standard quasiclassical approximation 



Wo = -Im 



K2 



am 



3^/3 



1TUJ 



2 K 2/3{z)dv, 



-Im 



2am 2 



3V3 



ituj 



3 



-i 



1 



v 



K 2/3 (z)dv, 



and the corresponding corrections to this approximation 

i 



W, 



(i) 



2 2 

am fi 



30V3 



TXUJK 



Gi(v,z)- 



dv 



(C.6) 



(C.7) 



where 

G 2 (v, z) = (36 + Av 2 - 18z 2 )K 1/3 {z) + {3v 2 - 57)zK 2/3 {z), 

G 3 (v, z) = -(34 + 2v 2 + 36z 2 )K 1/3 (z) + (78 - 6v 2 )zK 2/3 (z). (C.8) 

In derivation of Eqs. flC.7l) - flC.8j) the integration by parts was used as well as the 
intermediate equalities in Eq. (IC.5l) . The obtained representation of probabilities is 
convenient for calculation asymptotic at k ^> 1 which are 



w 3 



(i) 



2 2 

am jj, 



30^ 



-Wi 



w 2 



12Ak 1/3 - 90tt, 



TXUJK 



-11^ 1 / 3 + 84tt, A = 3^^^ = 8.191... 



(C.9) 



At k <C 1 the main contribution is given by the terms oc z 2 in Gi Eq.( 1C.8l) . One 
has 



lb 



(Mil 



! /i 2 2v/2 



2 ' luk 2 5V3 
The relative magnitude of the corrections at k ^> 1 



(i) 



2^ 2 

.,4/3 



2.9 



1/3 



(1) 



11/i 2 / 2.9 



9k 4 / 3 W/3 



(CIO) 



(C.ll) 
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At k <C 1 one has 

^- ~ (C 12) 

what agrees with Eq. (jB.7p . 



D Appendix 

On the derivation of created particle spectrum in SQA 

Substituting in Eq. (14.31) x = (1 + v)/2 one has 

dW SQA am 2 



dv ' 2V3 



1TLU 



2±±^tf 2/3 (*) + r K 1/3 (y)dy 

1 — V J z 



(D.l] 



where z is defined in Eq.f l3.2p . If one omits the integration over v in Eq. (l3.7p . the 
result of this manipulation disagrees with Eq.f lD.lh . So the described manipulation 
is erroneous. Only if one performs integration by parts of the second term in 
Eq. flD.lh and use the recursion relation 



zK 1/3 (z) = -z dJ ^ ~ §* a ,s(*). (D-2) 



one obtains Eq. (13.71) for integral probability. 
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Figure 1: The pair creation probability for unpolarized photon at /x = 0.1 as 
a function of r. The probability in saw-tooth form is calculated according to 
Eq. (12.1 21) . The dashed line is the result of averaging over the interval (— fi/2 +r -j- 
ji/2 + r). The solid line is calculated according to Eq. (13.41) . 
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2.5 5 7.5 10 12.5 15 17.5 20 

r 

Figure 2: The ratio of probabilities W {tK) /W {SQA) . The curves 2 and 3 are for 
polarized photons, the curve t is for unpolarized photon (see Eq. ( 13.41) and Eq. (13.71) ). 
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r 

Figure 3: The pair creation by a photon probability at \i — 3 in quasiclassical 
approximation. The upper curve is the probability in SQA Eq. (13.71) . The lower 
curve is the sum W^ SQA ^ + Eq.CT 
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5 10 15 20 25 30 



r 

Figure 4: The probability of pair creation at jj, — 3. The probability in saw- 
tooth form is calculated according to Eq. fl2.12p . The result of averaging over the 
interval (— /i/2 + r -4- fi/2 + r) is given by the dashed line. The solid curve, which 
is quasiclassical approximation, calculated according to Eqs. fl3.7p and fl3.ll) . 
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Figure 5: The spectrum of created particles for \i = 3, r = 110 (n m = 66) vs n (see 
Eq. fl4.2p ). The spectrum in saw-tooth form is calculated according to Eq. fl2.12l) . 
The dashed curve is the spectrum averaged over the interval (— 2/i + r -f- 2/i + r). 
The solid curve calculated according to Eq. fl4.3p . 
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Figure 6: The energy distribution of the created electron(positron) at /J, — 3 for 
different r: r = 30, n m = 16 (the upper curve), r = 110, n m = 66 (the middle 
curve), r = 500, n m = 318 (the lower curve), calculated according to Eq. (14.31) (the 
solid curve) and according to Eq. fl4.6p (the dashed curve). 
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